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ABSTRACT 

'  This  paper  la  concerned  with  the  analysis  of  the  solution  set  of  the  two- 
point  boundary  value  problem  modelling  the  avalanche  effect  in  semiconductor 
diodes  for  negative  applied  voltage*  -  This  effect  is  represented  by  a  large 
increase  of  the  absolute  value  of  the  current  starting  at  a  certain  reverse 
basis  We  interpret  the  avalanche-model  as  a  nonlinear  eigenvalue  problem 
(with  the  current  as  eigenparameter)  and  show  (using  a  priori  estimates  and  a 
well  known  theorem  on  the  structure  of  solution  sets  of  nonlinear  eigenvalue 
problems  for  compact  operators)  that  there  exists  an  unbounded  continuum  of 
solutions  Which  contains  a  solution  corresponding  to  every  negative  voltage* 
Therefore#  the  solution  branch  does  not  "break  down"  at  a  certain  threshold 
voltage  (as  expected  on  physical  grounds)*  We  discuss  the  current-voltage 
characteristic  and  prove  that  the  absolute  value  of  the  current  increases  at 
most  (and  at  least)  exponentially  in  the  avalanche  case  as  the  voltage 
decreases  to  minus  infinity* 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  we  investigate  the  mathematical  model  equations  for  impact 
ionisation  in  a  semiconductor  diode This  effect  (also  called  avalanche 
generation)  is  characterised  by  a  "sudden^  increase  of  the  current  flowing 
through  the  device  starting  at  a  certain  negative  voltage.  Physically,  the 
diode  "breaks  down"  shortly  after  the  onset  of  avalanche  generation. 
Therefore,  it  was  conjectured  that  there  is  a  threshold  voltage  beyond  which 
no  solutions  of  the  avalanche  model  exists.  We  show  that  this  conjecture  is 
false;  more  precisely  a  continuous  branch  of  solution  along  which  every 
negative  voltage  and  every  negative  bias  is  assumed  (at  least  once)  exists. 
Mathematically,  the  avalanche-effect  only  becomes  apparent  through  an 
exponential  increase  of  the  absolute  value  of  the  current  starting  at  a 
certain  negative  voltage. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sumary  lies  with  MR C,  and  not  with  the  author  of  this  report. 


A  NONLINEAR  EIGENVALUE  PROBLEM  MODELLING  THE 
AVALANCHE  EFFECT  IN  SEMICONDUCTOR  DIODES 


Ester  A.  Markovich* 

1.  INTRODUCTION 

No  investigate  tho  ( one-dimensional )  boundary -value  problem  which  describes  the 
performance  of  a  semiconductor  diode  in  the  case  of  avalanche  generation*  The  physical 
situation  is  as  follows*  A  semiconductor  is  doped  vlth  donor  atoms  on  the  right  side  In¬ 
side)  and  with  acceptor  atoms  on  the  left  side  (p-side)  and  a  bias  is  applied  to  the  Otaic 
contacts  (see  Figure  1)* 


pn- junction 


Anode  H 


r Cathode 


Figure  1 i  Diode 


For  simplicity  we  assume  that  the  pn-junctlon  is  in  the  middle  of  the  device,  that  the 
doping  profile  (that  is  the  difference  of  the  concentrations  of  donors  and  acceptors)  is 
constant  in  the  n-slde  as  well  as  in  the  p-slde  and  odd  about  the  pn- junction* 

A  well-known  phenomenon  is  the  'breakdown'  of  the  diode  due  to  impact  ionisation 
(avalanche  generation,  see  8se  (1981))  under  sufficiently  large  negative  bias*  This 
'breakdown*  is  based  on  a  'sudden'  increase  of  the  current  (as  a  function  of  the  applied 
bias)* 

To  study  the  current-voltage  (J-V)  characteristic  of  the  device  we  investigate  the 
basic  semiconductor  device  equations  describing  potential  and  carrier  distributions  In  the 


*Technische  UhiversitXt  Mien,  Znstitut  fur  Angevandte  und  Numerische  Ma theme t Ik, 
Gusshausstrasse  27-29,  a- 1040  Mien,  Austria,  Europe 


Sponsored  by  the  United  States  Army  under  contract  No*  DAAG29-80-C-0041 


diode  (see  Van  Roosbroeek  (1950),  Sse  (1981)); 
k2 


(1.1) 

(1.2) 

(1.3) 


X  8"  "  n  -  p  -  D  Poisson's  equation 

n'  •  n8*  ♦  alactron  continuity  equation 

p*  •  pf*  -  hole  continuity  equation 


•1  <  x  <  1 


8  denotes  the  electrostatic  potential,  8'  *■  the  electric  field,  n(p)  the  electron 
(hole)  density,  Jn(Jp)  the  electron  (hole)  current  density  and  D  the  doping  profile. 
The  equations  ( 1.1 )-( 1.3)  are  already  in  dimensionless  font,  the  doping  profile  is  scaled 
to  maximally  one  and  the  independent  variable  x  to  1-1, 1] .  In  our  symmetric  and 
piecewise  constant  case 


[  1,  0  <  x  <  1  (n-slde) 

(1.4)  0  *»  < 

(  -1,  -1  <  x  <  0  (p-side) 

2 

holds.  The  par  junction  is  at  x  ■  O.X  («  1)  is  a  scaling  parameter. 

Generally  the  current  relations  are  given  by 

(1.5)  (a)  J*  -  Ri  (b)  J#  -  -R 

n  p 

where  the  recombination-generation  term  R  is  a  nonlinear  function  of  n,  p,  Jn,  and 
8'  (the  electric  field).  We  assume  that  R  is  given  by  the  avalanche-generation  term 
(see  8se  (1981),  SchOt*  (1982))} 

(1.4)  R  -  R(J  .J  ,8* )  -  -«C8MC|J  I  ♦  |J  1) 

n  p  n  p 

where  a  >  0  is  the  electron-hole  ionisation  rate,  a  is  strongly  field-dependent. 

-  -2-  •  o 

Commonly  used  o*s  are  «(t)  -  ye  ^  ,  o(t)  •  title  Y,o  >  0.  FOr  simplicity  we 

assume  that  a  t  C((0,1])  ♦  (0,y)  ,  Y  >  0  is  the  nonnegative  functional 

(1.7)  o(f)  -  $(lfll-1  tJ),  8  i  (0,«)  ♦  (0,Yl *  $  €  C<  [<>,•})  and  nondecreasing 

(Ifl.  . .  i*  sup  If(x)l).  We  will  later  on  remark  on  the  extension  of  results  to  more 
[m,h]  a<x<b 
realistic  ionisation  rates. 

The  total  current  J  is  given  by 

(1.8)  J  • 

n  p 

Note  that  J  is  a  constant  in  (-1,1)  because  of  (1.5). 
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1h«  boundary  condition*  (at  tha  Ohmic  contacts)  tor  ( 1 . 1 )- (1 •  5)  ara 

(1.8)  (a)  np  -  n-p-D-0  at  x  <  i  1 

where  52(<  1)  also  originates  from  tha  scaling  and 


(1.8) (b) 


♦  d)  -  in 


n(  1 )  _ 


V,  *(-1) 


r 


V  e  A  is  tha  (scaled)  voltage  applied  to  tha  diode  (details  on  tha  scaling  can  be  found  in 
Markovich  and  Ringhofer  (1982)  and  Markovich  (1983)). 

Because  of  our  syesmtry  assumptions  ve  restrain  tha  investigation  to  'symmetric* 
solutions/  i.e.  solutions  vhich  fulfill 

(1.9)  *(x)  «  -  ♦(-*),  n(x)  -  p(-x),  Jn(x)  -  Jp( ~x) /  x  6  (-M)  . 

Another  simplification  is  accomplished  by  employing  the  substitution 

(1.10)  n  •  d2e*u,  p  -  «2e“*v,  J  -  d2e*u*,  J  -  -62e“*v'  . 

n  p 

The  system  of  equations  obtained  from  (1. 1)-(1.8)  by  using  (1.9),  (1.10)  is 

(1.11)  xV  -  «Vw  -  i2»~*v  -  1 

(1.12)  (.^u’)*  ■  -.(♦Mtl.^uM  +  |.  *v*|)  0  <  x  <  1 

(1.13)  (.  *v')’  -  -a(/')(|.*u'|  |.  V|> 

subject  to  the  boundary  conditions 


( 1* Id) (a) 


♦(0)  -  o,  0(1) 


,4!  J. IS) 

2tr 


( 1. 14)  (b) 


u(0)  •  v(0),  ud)  »  eV 


( 1. 14)  (c)  v'(0)  -  -u'10),  v(1)  -  e*V  . 

.y  V 

The  boundary  conditions  for  u  and  v  at  x  •  -1  are  u(-1)  •  e  r  v(-1)  -  e  .  The 
maximum  principle  (see  Protter  and  Weinberger  (1967)  applied  to  (1.12)/  (1.13)  gives 

(1.15)  u  3  e-lV‘,  v  >  .-<V|  on  1-1/1). 


Therefore  n  and  p  are  positive  (as  physically  required  for  densities).  A  solution 


for  V  -  0  la  9iv.11  by  u  5  1,  v  1  1  and  b y  solving 


rr- 


(1. 16)(a) 


( 1. 16) (b) 


- « V  -  «V*  -  i,  o  <  x  <  i 


♦  (0)  -  0,  #(1)  -  Inf1-—  1  *  ) 

262 


th®  solution  (V,*,u,v)  -  (0,*^,1,1)  where  ♦  is  ths  unique  solution  of  (1.16))  is 
csllsd  equilibria  solution.  It  isplies  J  -  0  (the  whole  diode  is  in  thermal 
equilibrium). 

*>•  two-point  boundary-value  problem  (1. 11)-(1. 14)  models  the  bias-controlled  diode. 
In  some  cases  it  is  more  convenient  to  investigate  the  current-controlled  device 
represented  by  the  equations  (1. 11 )-(1. 13)  subject  to  the  boundary  conditions. 


( 1. 17) (a) 


( 1. 17) (b) 


(1.17(c) 


♦  (0)  -  0,  +  (1)  -  inf1  *  ■  -1/  )  -  in  u(1) 

2r 


u(0 )  -  v(0),  u(1)v(1)  -  1,  u(1)  >  0 


u'(O)  -  v'  (0 )  -  -  -~ 

26*  26* 


(note  that  (1.17)(c)  follows  from 

J  -  J  (0)  ♦  J  (0)  -  62.*(0)u*(0)  -  62."*<0)V'(0)  -  262u’(0)  -  -262V*(0)>  . 
n  p 

The  problems  (1.11)- (1.14)  and  (1.11)- (1.13),  (1.17)  are  equivalent  in  the  following 


sense.  A  solution 


(Vj #Uj ,v^ )  of  ( 1. 11 )-( 1. 14)  yields  the  solution  (J1 ,u^ ,v^ ) 


( 1. 11 )-( 1. 13),  (1.17)  where  *  62e  *uj  -  62e  *v*  and  a  solution  ^2f^2,U2,V2^  °* 

(1.11)- (1.13),  (1.17)  yields  the  solution  <v2'*2'u2'v2)  of  (1.11  )-(1.14)  where 

V2  "  tn  tt2(1)* 

There  are  numerous  analytical  and  numerical  ivestigations  of  the  (even  multi¬ 
dimensional)  semiconductor  device  equations  in  the  non-avalanche  case  (i.e.  the 
recombination-generation  rate  R  only  depends  on  n  and  p)  (see  Nock  (1983)  for  a  rather 


r 


I 


complete  presentation  of  the  results  es  well  as  for  e  collection  of  references)*  For  the 
avalanche  problem  however  there  are  (to  the  author* a  knowledge)  only  a  few  numerical 
studies  (see  SchOts  (1982);  Schtttx,  Selberherr  and  ftStsl  (1962)), 

In  this  paper  we  regard  ( 1. 1 1 )-( 1, 14)  and  (1.11 )-(1. 13),  (1.17)  as  nonlinear 
eigenvalue  problems  (in  the  eense  of  Kabinowlts  (1971),  Krasnoeelskii  (1964)  and 
investigate  the  solution  set  for  nonpositive  current 

n.ie>  c"  -  {«,♦,,,»)  e  (—,o]  *  (c2< [0,i] )3| (♦,.,▼) 

•olv..  (1. 11)-(1 . 13),  (1.17)  with  J  -  I) 

and  the  properties  of  the  cur rent- voltage  (J  -  V)  characteristic 

(1.19)  *  {(v,j)  e  ft  *  (-•♦Ol)  there  is  (4,u,v) 

such  that  (J,4,u,v)  e  c  and  V  •  in  u<1)}  . 

the  wain  theoreai  of  this  paper  states  that  c"  contains  an  unbounded  continuu»  (i.e.  a 

2  3 

closed  and  connected  set  in  the  (-*•,<>}  *  (C  ([0,1]))  -topology)  emanating  fro*  the 
equilibrium  solution  (0,4^,  1,1)  whose  projection  into  (-<»,0)  equals  (-*#0]  (that 
means  C~  contains  solutions  for  all  J  <  0)  and  that  the  voltage  V  ♦  as  J  ♦ 
Therefore  ( 1, 11 )-( 1. 14)  has  a  solution  (4#u,v)  for  every  V  <  0. 

This  result  holds  independently  of  the  upper  bound  y  of  the  ionization  rate  a  and 
carries  over  to  more  realistic  ct's  than  given  by  (1.7).  Therefore  the  conjecture  that 
the  branch  of  solutions  of  (1.11)-(1. 14)  breaks  down  if  Y  >  J  holds  (see  See  (1981))  is 
mathematically  rejected  at  least  for  this  model  problem.  We  show,  however,  that  the 
magnitude  of  y  has  a  decisive  impact  on  the  J  -  V -characteristic.  For  a  5  0 
(nonavalanche  case)  the  current  fulfills  c^V  <  J  4  c^V  for  V  <  0  while  |J|  increases 
exponentially  as  V  ♦  -•  for  Y  >  J  (c1#c2  >  0  only  depend  on  X  and  «>.  The 
exponential  growth  of  the  current  represents  the  * avalanche  effect'  and  the  diode  'breaks 
down'  in  real  life  when  the  current  gets  too  large.  We  also  show  nonuniqueness  for  V  *  0 
for  all  a  for  which  a(4^)  is  sufficiently  large* 

The  paper  is  organised  as  follows.  Section  2  deals  with  the  a  priori  estimates  needed 
to  prove  existence  of  solutions  for  all  J  <  0  and  Section  3  contains  the  existence  proof 


and  conclusions 


2.  A  PKIORI  KSTIKATM 


For  ths  following  wo  taka  J  <  0. 

At  first  ws  solvs  ths  continuity  aquations  (1 • 12)-( 1. 13)  for  flxsd  f  C 
rswrlts  thssi  as 


<2. 1 ) (a) 


(2. 1 )  (b) 


J*  -  -o(t')(|J  |  ♦  |J  |) 
n  n  p 


J'  -  <*(♦'  )(|J  I  ♦  |J  I) 
P  «  P 


0  <  x  <  1 


with  ths  initial  valusa 

(2.2) (c)  J  (0)  -  |  ,  J  <0 )  “  “  • 

n  /  p  c 

(2.1) (a)  impllss  that  is  nonincraasing*  sines  JR(0)  <  0  ws  gst  <  0  on  (0,1 J, 

Jp<0)  <  0  holds  and  thsrsfors  ws  (initially)  solvs 


J*  ■  -a|  J| ,  J*  «  a| J| 
n  p 


(ws  nostly  drop  ths  argument  ♦  '  of  a)  and  gst 


(2.3) (a)  J 


•  -  (2ox  +  1),  J  -  -  <-2ax  ♦  1)  for 

2  p  2 


x  e  [0,1]  if  0  <  a  <  j 


*  *  lt  *  >  j 


For  a  >  “  ws  havs  to  solvs 


J*  -  a(J  -  J  ),  J1  -  -o(J  -  J),  —  <  x  <  1 
n  n  p  p  n  p  2a 


j  -  -|j|,  j  (^-)  -  o 

n,'2o''  pv2a' 


and  obtain 


(2. 3) (b)  J  m  m  (s2^1  ♦  1),  J  -  -  ^  (-s 
n  2  p  2 

u  and  v  ars  computed  from  (1.10),  (1.14){b),(c)s 


2ax-1 


1)  for  x  0  l]  . 


-6- 


I 


(2.4)<a)  u  ■  ®V  +  • 

2«2 


/  •”****  (2c»  +  t)d®  for  x  e  [0,1]  If  0  <  a  < 

v  < 


^^  (•) ,  2<*g-1 


1/2a 


'<• -  + 


1)d«  for  x  e  [^,  l]  if  a  >  - 


(2a®  +  1)ds  +  /  ®"*<#)<®2o*“1  ♦  l)d* 


1 

1/2o 

for  x  «  [0,  j]  if  a  >  -1 


r  i  .(  . 

I  +  1)d»  for  x  e  (0,1]  If  0  <  •  <  | 

X  2 


,-v  .  JJL 
2«2 


♦  (®),  2a®-1 


♦  1)ds 


for  x  e  l]  It  a  >| 


j,/2“  _2a.  +  1)d.  +  j1  +  nd> 

x  1/2a 


for  x  «  [«,  j]  if  a  >  | 


W®  us®  th®  condition  u(0)  -  v(0)  to  r®l«t®  V  and  J  and  g®t 


(2«  5) 


®inh(J*I^>),  J  <  0 
4«Z 


where  th®  functional  I  :  C  ([0,1])  ♦  R  i»  given  by 


(2.6)(a) 

with 


!<♦>  -  /  (••‘“’g  (■>  ♦  „(•)!«• 


(2.6)(b)  g  (x)  - 
<s 


-2ax+1  for  x  e  [0,1]  if  0  <  a  <  \  and  for  x  e  [o,  —]  if 


2ajr*1 


1  for  x  *  11  if  «  >  4 


-7- 


I 


r 


(  2ox+1  for  x  €  [0,1]  if  0  <  a  <  1  and  for  x  e  [o,  1^]  if  a  > 


(2.6)<b)  fa(x)  - 


2ax-1 


♦  1  for  x  6  [i^,  l]  if  a  >  ^ 


I  e  0(0* ( [0,1] )  ♦  R)  holds.  For  the  estimates  of  the  current  J  in  terns  of  the 
voltage  V  we  use 
(2.7) 


a  .  1  lf  1<t)  *  o 


!(♦) 

We  collect  the  properties  of  u,  v,  Jn,  Jp  and  J  in: 

2.1 :  Assume  that  0  <  a(6'  )  <  ^  holds.  Then 
1  *  1 

(i)  1(f)  >  (1  -  2a)  /  e*t8)ds  +  /  e*"^^*)ds  >  0 

0  0 

(ii)  J  <  0  <—  >  V  <  0,  J  -  0  <— >  V  -  0  <— >  u  =  v  =  1. 


(ill)  Let  V  <  0  hold. 

Then  J„  <  0,  Jp  <  0  on  (0,1) 

(iv)  u  is  decreasing  on  [0,1],  eV  <  u  <  e  V 
(v)  v  is  Increasing  on  (0,1],  eV  <  v  <  e  V. 
2. 2:  Assume  that  ®(^* )  >  ^  is  fixed.  Then 


(i)  there  is  a  ♦  f  c\[0,1]>  such  that  1(f)  <  0 


(ii) 

O 

0 

> 

A 

1 

o 

1 

*1 

<—> 

e 

<ii 

< 

hi 

i  1),  J  <  0  <—  >  (V  < 

0,  !(♦)  >  0  or 

(iii) 

V  >  0,  I(t)  <  0) 

V  •  0  <— >  J  »  0 

or 

!(♦)  - 

0. 

(iv) 

Let  J  <  0  hold. 

Then  Jn  <  0  on 

[0, 

11;  Jp 

<  0  on  [o,  ij),  Jp 

10 

c 

0 

o 

A 

(v) 

u  is  decreasing 

on 

(0,1], 

u  >  eV 

(Vi) 

v  is  Increasing 

on 

1°'  sl 

and  decreasing  on  [ 

h-  1].  *>*v- 

m 

0 

* 

A 

for  any  solution 

(V  - 

0,**,u* 

»v* )  *  (V  -  0,6  ,1,1) 

of  (1.11)- (1.14) 

Y  >  j)  !(♦*)  -  0  has  to  hold. 

We  now  turn  to  Poisson's  equation  (1.11)  subject  to  the  boundary  conditions 
(1.14) (a).  Differentiating  (1.11)  gives 


(2.8 )  (a) 


>"  -  (n  +  p )f  +  J,  0  <  x  <  1 


and  (1.14)  implies 


(2.8) (b)  »  -  -L,  (♦’(1))'  -  0  . 

xz 

The  maximum  principle  (see  Protter  and  Weinberger  (1967))  yields 

Lemma  2.3:  J  <  0  Implies  f  >  0  on  [0,1]  and  V  <  f  (1),  Thus  the  a  priori  bounds 

(2.9)  0  <  *(x)  <  *<1)  -  *e(D  -  v,  X  e  [0,13 
follow  from  J  <  0. 

Differentiating  (2.8) (a)  and  using  (1.11)(a)  yields 


( 2. 10  )  ( a) 

A2(*")"  -  t(n 

♦  P>  ♦  *V2J**  ♦ 

{*• >2  +  (J  -  J  )♦’,  0  <  x  <  1 

n  p 

( 2. 10 ) (b) 

*"<0)  -  -  L, 

JT 

♦"(1)  -  0  . 

From  (2.3) 

we  conclude  that 

Jn  -  Jp  <  0  in 

[0,1)  for  J  <  0.  Since  >  0  we  obtain 

Lemma  2.4. 

J  <  0  implies 

♦"  >  -L  on  [0,1] 
A 

and  therefore 

(2.11 ) 


n(  x)  >  p(  x) ,  0  <  x  <  1 


holds. 


Proofs  *  *  -  ~  is  a  lower  solution  of  (2.10)  and  (2.11)  follows  from  (1.1)  (with 
~ —  X2 

D  =  1  on  [0,1] ). 

We  now  derive  upper  bounds  for  n. 

Lemma  2.5.  Let  J  <  0  hold.  Then 


(2.12, 


{2a  +  1 ) ( 1  -  x )  for  x  e  [0,1]  if  0  <  a  <  j 
(e2a_1  +  1)(1  -  x)  for  x  e  [Lj,  l]  if  a  >  j 


2 

2 


2x  ♦  <e2a-1+  1)(l  -  ^) 

for  x  e  [C,  jj]  if  a  >  -j  • 


Proof.  We  multiply  (2.4) (a)  by 
s  >  x  ( f  is  nondecreasing). 

In  the  case  of  zero  generation 
independent  of  J  and  V. 


62e^  (getting  n)  and  estimate  <  1  for 

(2.12)  is  then  obtained  by  integration  and  estimating 

—  (1  -  .-2a<1-x)>  <  1  -  x  . 

2a 

(a  =  0)  we  obtain  upper  bounds  for  n  and  p  which  are 
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La—  2.6.  Lat  a  3  0  and  J  <  0  hold.  Than  *■  <  0  on  [0,1]  and 

(2.13)  P(x)  <  n(x)  <  p(x)  ♦  1,  0  <  x  <  1 

(2.14)  0  <  n(x)  +  p(x)  <  /  1  +  4«4,  0  <  x  <  1 

hold. 

J 

Proof .  o*0  implies  J  5  J  =  —  and  therefore  z  *  0  is  an  upper  solution  of 

“  1 "™  n  p  2 

(2.10).  (2.13)  follows  from  (1.1)  (a).  Also  (n  +  p) '  -  ♦  '  (n  -  p)  holds*  Thus  n  ♦  p 

nondecreasing  and  n  +  p<n(1)  +  p(U-  /7T  46*  follows. 

We  now  derive  lower  and  upper  bounds  for  +  using  the  estimates  of  n  and  p. 
Tafia  2.7.  Let  J  <  0  hold.  Then 

-  2 

(2.15)  fix)  <  [f(1)  ♦  —  )x  -  *-7,  0  <  x  <  1 

2A2  2A 

follows.  If  0=0 

2 

(2.16)  *(x)  >  -  ♦  ♦(l)x 

21Z 

holds • 

1  2 

Proof.  i(x)  -  (♦<!)  ♦  —7 )x  -  *-r  solves 
2A2  2A 

12*«  -  -1,  *<0)  -  0,  *(1)  -  *(1)  . 

leona  2.4  implies  that  f  is  an  upper  solution  of  (1.11)  (a).  For  a  =  0  we  obtain  by 

integrating  ♦"  <  0  :  pi 1 )  <  (0 )  and  from  ~ 

A2 

2  2 

*(X)  >  -  ♦  ♦' (0)x  >  -  *-t  ♦  *(i)x  . 

212  21Z 


To  get  a  lower  bound  for  p  in  the  avalanche  case  we  prove 

Lena  2.8.  Let  j  <  0,  |J|  >  K1  hold.  Also  let  0  <  e  <  1  be  such  that  |j|e5  >  K2 
where  K^#  Kj  only  depend  on  A,  6  and  Y*  Then 


(2.17) 

holds  where 


♦  •(x>  >  .  ■1  T7  *  ®(*>  1 - ,  0  <  X  <  1 


e  ♦  ga(x) 


IJI* 


lol 


[01]  *  *3  ^K3  only  depends  on  A ,6  and  Y)  and 
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is 


(2.18) 


(2<*  ♦  1)(1  -  x)  for  x  e  (0.1)  If  0  «  •  <  j 


g0(x) 


(«2*  1  ♦  1  > (1  -  x)  for  x  6  [|j,  l]  if  «  >  j 
j  -  2x  >  I.2®"’  ♦  1)(1  -  for  x  8  [0.  if  «  >  i  . 


Proof.  The  Tiwm  2.4  end  2.5  imply 


2 S_ 

Ml 


< 


♦  Z77 


44* 


Ml 


♦  9a(*> 


We  now  choose  €  euch  that  1  ♦  /  1  ♦  44*  <  |J|c.  Thu*  <  c  ♦  90(*>  holds. 

Obviously  the  solution  y<>  0)  of 


(2. 19) (e) 


(c  ♦  g  (x))y  -  1 
o 


(2. 19) (b) 


y#(0) 


y'd) 


o 


is  a  lower  solution  of  (2.8),  that  means  0  <  y  <  9*  holds  on  [0,1].  For  large  |J|  the 

problem  constitutes  a  linear  singularly  perturbed  (Neumann-'type)  boundary  value  problem 

X2 

with  the  reduced  solution  (obtained  by  setting  jyp  to  sero) 


1 

yr  *  e  ♦  g  (x)  * 

Cl 

A  standard  singular  perturbation  analysis  (see  Howes  (1978))  which  takes  the  possible 
smallness  of  t  into  account  <y ^ ( 1 )  *  “  II)  gives 


,y  ■  yr* ro, 11  *  *3 


/  |J|«* 


whenever  |J|  >  Kj,  |j|c  >  This  implies  (2.17). 

A  lower  bound  for  +  follows  by  integrating  (2. 17 )i 
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2.9.  Let  tha  assumption  of  Leama  2.7  hold*  Than 

J1  *  STTT  An  ♦  (2a  ♦IMI  -  x>) 


(2*20) 


9(x)  >  Mix) 


|J|« 


for  x  e  (0,1)  if  0  <  a  <  - 
l2  for  X  e  [o,  n  a  >  J 


L,  ♦  — — - 
3  2a-1 


*n  (- 


2a*1  * 

♦  1  (a  ♦  D(1  -  x)  ♦  e 


for  x  e  [*jj,  l]  If  a  >  i 


holds  whsre  1)  <  *3  M<1  LV  **2#  L3  onlY  dspsnd  on  Y- 

Sines  c  can  ba  aads  arbitrarily  mull  whan  J  ♦  (still  kaaping  |J|c5  large) 
th*  Uaui  2.8  and  2.9  ixply  that  t'lD  >  *4  and  ♦<1>(>  Kj  tn  j)  b*coaw 

unboundad  as  J  ♦  -*•- 

Ha  also  naad  an  uppar  bound  of  ♦’ : 

Lass  2*10.  J  <  0  implies 

(2.21)  (0  <>♦'  <  K6*|V,(|J|  ♦  1) 

vtiara  only  dapands  on  X  and  6* 

2  V 

Proof i  (2. 4) (a)  lapllas  n+p>n>6a.  Thus  tha  solution  w  of 

X2w"  -  «2*Vw  ~  |J|,  *’(0)  -  -  »'(D  -  0 

%  * 


is  an  uppar  solution  of  (2*8).  There fora 


4f(x)  <  w(x) 


«2 


cosh  (j  aV^(1  -  x)) 
X«aV/2.inh(f  aV/2) 


and  (2.21)  follows  sinca  v<  4^(1). 

Na  now  employ  tha  derived  bounds  to  gat  a  priori  estiaates  on  tha  current-voltage 
characteristic. 

Thaoraa  2. 1.  Assume  that  <J,t#u,v)  e  c  end  (V, J)  e  J~.  Than 
(2,22)  Cjvl  <  |J|  <  C2|V|  if  a (9*)  -  0 
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(2.23) 


c,|v|  «  |j|  <  c3«|v|  (1  -  •_J|V|)  If  0  <  «(♦•)  <  | 


holds.  If  «<♦* )  >  -  thon 


(2.24) 


c4<.  -  .-i|,S."ln  •  <  |J|  <  fi  „.,,s  . 


holds  for  j  <  «(♦')  <  aQ  with  snw  oQ  <  1  snd  every  U  >  0.  If  a(^')  >  thon 

(2.25)  Cjd  -  ."2,V,).,V|  <  |J|  <  D  «xp<<5  ♦  MX.2*'1  ♦  1)|V|)  . 

Ths  constants  C1#...#CS  only  dapsnd  on  X#6  and  y.  0  depends  on  A,5#y  and  p. 

Proof.  0  <  a  <  -j  and  ♦  >  0  imply  0  <  1(f)  <  2  /  cosh  f(s)ds  and  (2. IS),  (2.7)  ylaldt 
4  0 

2  alnh  (♦<!)  ♦  -^r) 

*«  • 


0  <  1(f)  4 


A hi 


♦CD  ♦ 


2A 


Thus 


«2<i  -  .“2|v,)(*(i>  ♦  -*r) 


ui  > 


.',V,.li»h(t(1)  ♦  ^r) 

2Xa 


lil 


and  ths  lower  bounds  for  |J|  in  (2.22),  <2. 23)  follow,  a  t  0  implies 
I 

I (♦)  ■  2  /  cosh  ♦(s)ds  and  the  estimate  (2.16)  glees 
0 

•l»h(+(1)  -  -ij) 

!<♦)  >  2  - 

♦O)  - 

2Xa 

when  |V|  is  so  large  that  fd)  *  ~~  >  0.  We  derive  (using  (2.7)) 

2XJ 


UI  < 


«2(1  -  .'2,TI  )(t(1)  -  -lj) 
_ _ 2\ 

."|V,.l«h(*(1)  - 

2X* 


and  (2.22)  is  proven.  For  0  <  «  <  j  we  estimate 
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Thus 


!(♦>  >  /  (-2m  ♦  DM  •  1  -  a  . 

0 


m2|V|  -2|V|. 

Ul  <  *  *  \\  * - L 


and  tha  uppsr  bound  in  (2.23)  follows*  Mow  lat  a(4#)  >  -.  Than,  alnoa  4  la 
nondacraaaing  and  alnca  gft(x)  1»  positive  in  [o,  —)  and  negative  In  (j,  l]  i 


1(4)  <  /  f  <s)ds  +  f(1/2a)  /  9  (s)ds 
0  0  a 


4o  2a 


<fa  >0  in  CO#  tJ  f)  holds.  Ths  function  h(a)  -  1  *  4a  ~  2a  ***  ^  *  aniqus  ssro 


*0  e  ^2#  Thu* 


1(4)  < 


2  ,  a  >  «„ 


«  1^1  2a-1  A  1  +(1/2a)  1  ^ 

1  -  — —  ♦  ~  a  ♦Ta  $  —  <  a  <  a 

4a  2a  2  2  0 


holds*  Zf  1(4)  >  0  ths  lower  bound  in  (2.25)  follows  and  ths  lowar  bound  in  (2*24)  (also 
for  1(4)  >  0)  Is  iapliad  by  (2.15)  which  gives 

♦&) «  c  *  |v|  h  • 

lass  2*9  inplias  that  'V  ♦  as  J  ♦  and  tharafora  1(4)  <  0  can  (for  J  <  0) 

only  hold  for  |J|  <  F,  whara  F  dapsnds  on  1,4  and  y.  This  proves  tha  lowar  bounds 
in  (2*24),  (2.25). 

(2.20)  yi.lda  fO)  •  *  (11  -  V  >  t  ♦  -r-| -  tn  7.  M  Mt  «  -  1/1  J| 1/5"®  for 

a  j  .  t  c 

1  •  ♦ 1 

0  <  o  <  —  and  obtain  th.  upper  bound.  In  (2.24),  (2.25)  for  |J|  auff.  larq.  alne. 


V  <  ♦  (1)  holda. 


Sines 


2.S  inplias  that  a(4*)  ♦  y  as  J  ♦  -•  (or  as  V  ♦  -•).  Tha  thaorsa 


prows  that  tha  currant  J  incraasas  (in  absoluta  value)  at  laast  (and  at  nost) 
exponentially  as  V  ♦  -•  in  tha  avalanche  casa  Y  >  For  saro  generation  (y  •  0)  th# 
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incraaa.  l«  at  aoat  (and  at  laaat)  llnaar.  For  tha  intarmadlata  eaaa  0  <  y  <  j  tha 
lncraaaa  la  at  laaat  llnaar  and  at  no at  axponontlal.  Tha  author  oonjaeturaa  that  tha 
diatlnctlon  of  tha  eaaaa  J  <  Y  <  o„  and  y  >  aQ  only  coaos  In  for  taehnloal  raaaona  and 
that  (2.25)  holda  for  all  Y  > 


3.  MISTKWC1  THKMBMS 

We  need  the  following 

Low—  3.  U  let  v  bo  tho  (unique)  solution  of  the  problow 
<3.1)(«)  w"  -  o(x)oW+n  -  b(x)d^f"1*  ♦  f(x),  0  <  x  <  1 

(3a1)(b)  w(0)  -  y^,  w(1)  •  y^ 

whoro  a,b  €  C(  (0, 1)  )r  a,b  >  0  on  (0,1)  f  Ml.  Then  w  »  w(n,y0#y1  ,f  ,•  )  regarded  oo 

napping  froa  *3  *  C1  ( (0#  1] )  into  C1<(0,1])  ia  completely  continuous, 
ftroofi  we  take  n  «  (3#nl,  yQ  e  [|«0  #y q)  #  y^  €  [\1 1  r  y  1 )  and  denote 

a  win  a(x),  a  i-  wax  a(x)  (analogously  for  b).  Then  the  unique  solution  w,  of 

0<x<  1  0<x<  1  1 

-  *^r!i 

*? «< *< 1 
*i<0>  “  Ho'  wi(1>  ”  !*! 


la  a  lowar  aolutlon  of  (3.1)  and  tho  uniqua  aolutlon  «2  of 

,»  0  <  x  <  1 

IV,  I  I 

*2(0)  *  V  V1*  *  “l 


_  -w.-n 

*2  "  ”  ,f*  10,1] 1 


la  an  uppar  aolutlon  of  (3.1),  l.a.  on  [0,1]  holda.  dlnoa 


.  -  "»i“3 

'•’•to.i)"-  ♦»*  ♦,f,(0.1] 

holds,  Ascolivs  theorew  implies  that  w  i  *  C*((0, 1))  ♦  c\[0,1])  nape  bounded  sets 
Into  precowpact  sets.  The  continuity  of  w  is  lwwsdiate. 

Now  we  prove  the  basic 
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Th ooron  3Jt  for  «ny  Y  >  0  tho  solution  got  C  of  (1. 11)-(1*  13),  (1.17)  oontalns  on 

■ting  from  tho 


~~  2  3 

unbounded  continuum  C  Un  tho  (^#0]  *  (C  ((0,1)))  -topology) 


equilibrium  ooltttlon  (0#9^#1,1)  whose  projection  into  (-«,0]  equals  (^*,0)  (l.o.  C 


contolno  o  solution  (J,#,u,v)  for  every  J  <  0). 

Froofi  W.  r«9«rd  V  -  V(J,#>  (firm  by  (3.5>>  a.  functional  V  i  (-*,0]  a  c’jKMJ)  ♦  m. 
Tho  continuity  of  Z  loplloo  tho  continuity  of  V.  doing  (2.3)- (2.5)  wo  rowrito  tolseon's 
oquntlon  (1.11)  oo 

(3.2) (a)  X2*-  -  .  g3m~+~VlJ'*)  _  |  +  |J|G(*)(x>,  0  <  x  «  1 


(3.2)(b) 

with 


♦  (0)  -  Of  0(1)  -  0(1)  -  V(J,#) 


(3.3) 


.<»><»  -  /’  |J^(  «.  .  /'  .*«•>  |i£|  « 


(x).  (x). 


Ini  Id  I 

“j  “I,  I  j  |  oro  indopondont  of  J#  thoy  only  depend  on  «(♦*)  and  on 
«>•  6  t  C1((0#1])  ♦  C1((0#1])  in  continuous  slnco  o  t  c\(0,1])  ♦  (0,y1  lx  continuous 

«o  got  t  -  ♦  ♦  4  ond  rowrito  (3.2)  no  tho  flood  point  probloo  +  •  T(J,9)  whoro 


y  -  T(J,«)  is  doflnod  so  tho  unique  golution  of  tho  probloo 


(3.4) (a)  X2y«  -  <2«xp(y  ♦  «  +  V(J,f  ♦  .)) 


-  «2  Mfl-j  "  ‘  »M,»#  ♦  a)>  -  1  -  X2*"  ♦  |J|0(f#  ♦  a),  0  <  x  <  1 


(3.41(b) 


y(0)  -  0,  y(1)  -  -*(J,*  ♦  a)  . 


3.1  ond  tho  continuity  of  V  and  0  imply  that  T  t  (—,0)  *  c  ((0,1))  is 
lotoly  continuouo.  v(0,w)  l  0  ond  thoroforo  y  "  T(0fc)  is  gluon  by  tho  golution  of 


4  )  J  9  "7"^.  4 

\V  -  «  •  -  «  o  •  -  1  -  X  0  <  x  <  1 


(3.5) (o) 

( 3.5)  (b)  y(0)  -  y<  1 )  -  0  . 

y  l  0  follows,  freo  Hob inowl tt  (1971,  Thooroo  3.2)  wo  conclude  that  tho  golution  sot  of 
T(J,g)  •  9  oontalna  an  unbounded  continuum  t”  (in  tho  (-«,0]  *  CY(  [0, 1])-topology)  • 
Thoorom  2.1  and  Lomno  2.3  imply 
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|V|  <  ■•*(♦  n>, 

•  C1 


and 


\0)  -  \  *  ivl 


2.10  yields 


♦;<♦•<  v 


D  -  ♦: 


ifa  ooncloda  froa  thaaa  eatiaatee 


Ul/c  Ul/c, 

10,11  *  10,1]  *  "t1  ♦  1,1  *  •  1*1) 


ahere  N  >  0  la  Indapandaat  of  J  <and  V).  tinea  S~  la  aaboandad  It  hat  to  ooatala 
aolutlone  (J,4)  for  all  J  <  0.  The  atateaent  of  tha  thaoraa  folloaa  by  ofata  rainy  that 
a  and  a  at  given  by  (2.4)  ara  ooatinaoua  at  faaetloat  of  (J##)  In  tha 
<-•,0]  *  C1 1  f0#  11 )  -topology* 

Tha  aoat  important  iapllcatlon  of  Thaoraa  2*1  lai 
Thaoraa  3.2i  For  any  T  >  0  tha  curran  t-vol tag*  char  act  aria  tic  J*  oontalna  a  oootinuoue 
curva  r"  aaanatlng  froa  (0,0)  ahoee  projection  Into  tha  J-axle  aqua  la  (^,0)  and 

ahoee  projection  I*“  into  tha  9-axle  fulfill* 

(3*6)  (a)  -  <-*,0)  for  0(y  <J 

( 3.6) (b)  (— ,0]  £fjs  (— for  Y  >  j  * 

Proof  i  Thaoraa  3*1  lapliaa  that  oontalna  a  oontlnuua  O*  aaana ting  froa  (0,0) 

ahoaa  projection  Into  tha  J-axla  equal  a  (  —  ,0].  Proa  Laaaa  2.9  an  ooncloda  that 
+  (1)  -  -  g  i»  poeltlve  and  unbounded  at  *  ♦  — .  Xaaaa  2.1,  (11)  lapliaa  that 

V  <  0  for  J  <  0  and  0  <  Y  <  ^*  Tbhrefore  (3.6)(a)  folloaa*  (3.6)(b)  la  concluded  by 
noting  that  *(1!  *  *#(1>  -  9  >  0  for  J  «  0  holde  (tea  laaaa  2.3). 

Tha  obvloue  conteqaenoe  for  tha  eolation  net 

(3.7)  o’  -  {(u,t,u,v)  «  (-*,0)  »  (C2((0,lj))3|(9,u,y) 
eolvee  (1. 11)-(1.14)  for  9*0) 

of  the  voltage-controlled  diode  la 
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Corollary  3.1  d”  contains  an  unbounded  oontlnuus  o~  containing  (0,t*,v*,u*»  whose 

projection  Into  (—.0]  equals  <—,0)* (♦•,«* ,y*)  -  (*^,1,1)  if  thn  equilibria  solution 

1  .. 

is  unique  (e.g*  for  0  <  t  <  Jl.  J  «  0  holds  for  every  (V,p,u,v)  •  o  .  (0  contains  a 

solution  for  every  v  <  0). 

W»  now  show  that  wultiple  solutions  of  (1.11)-(1.14)  for  V  -  0  can  occur. 

™»sor—  3.3i  Assusw  that  o(pM  >  ■(»  j)  where  H  >  0  only  depends  on  1  and  «.  Then 
there  is  a  solution  (♦*,u»,v»)  of  (1.11)-(1.14)  for  V  -  0  which  is  different  fres  the 
squilibriun  solution  1,1)  and  J*  ■  42s^*(u*)*  -  4Js”**(v»)'  <  0  holds. 

Proof i  wo  obtain  from  (2.6) 

t  1 

I(t)  ■  2  (/  cosh  +<s)ds  -  J  h  (•)  ainh  f<s)ds) 

0  0  a 

whsrs 


ha(x>  -  i 


2ax,  0  <  x  < 


i_  <  x  <  , 

2a 


(for  a(f)  >  . 


Obviously  h  (x)  >  2ax  on  [0,1]  and  tharafora 
<s 

1  1 

!<♦>  <  2  (/  cosh  +(s)ds  -2a  /  s  sinh  t<s)ds) 
0  0 


holds.  Choosing  a  such  that 


1 

/  cosh  f  (s)ds 
- : — 

/  a  sinh  f  (s) ds 

0  • 


iapliaa  «♦  >  <  0* 
I<*#  ♦  ♦)  <0  for 
Sines  tha  continuum 
intarsaction  of  E~ 
Of  (1.11)-<1e13),  (t 


Thus  thara  is  a  naighbourhood  W  of  0  in  cNlOpI])  with 
♦  ON.  (2s S)  implias  ♦  ♦)  >  0  for  9  [-»,0)  x  w. 

(usad  in  tha  proof  of  Thaoram  3s 1)  smanatss  in  (0,0 )•  Tha 
and  (-»,0)  x  h  is  not  smpty.  This  iapliti  that  thara  ara  solutions 
17)  for  J  <  0  for  which  V  >  0  holds*  Sines  V  is  nsgativa  for 
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J  <  Ot  |J|  sufficiently  large,  there  hee  to  be  6  *  t  J  *  0,  with 

♦  ♦*)  "  0*  This  gives  +  4*)  •  0  end  Theorem  3.3  follows. 

The  condition  Z(^)  <  0  implies  thet  J“  is  not  contained  in  (•*•,(>)  ^  (see 
Figure  2). 

J 


V 

Figure  2.  Qualitative  structure  of  the  J  -  V-characterlstlc  for  I(t#)  <  0. 

However,  I(^#)  <  0  is  physically  unreasonable  and  it  is  not  clear  whether  the 
nonuniqueness  of  the  equilibrium  solution  prevails  if  I(f^)  >  0. 

We  conclude  from  the  Theorems  2.1  and  3.2  that  the  avalanche  case  y  >  y  is 
distinguished  from  the  nonavalanche  case  0  <  y  <  -j  by  a  more  rapid  decrease  of  the 
current  J  as  V  ♦  -•  (see  Figure  3). 


Figure  3.  Qsalltative  Structure  of  the  J  -  V-characterlstic  for 
various  y's  ( I C )  >  0  is  assumed)  for  reverse  bias. 


Me  remark  that  an  investigation  of  (1*11)- (1*13)/  (1*17)  for  nonnegative  currant  can  be 
dona  in  a  similar  fashion*  A  ralation  of  tha  term  (2,5)  (with  a  diffarant  functional  1) 
holds  and  tha  axistanca  of  an  unbounded  continuum  of  solutions  E+  follows  as  in  tha  proof 
of  Theorem  3.1*  To  conclude  that  V*  contains  solutions  for  all  J  >  0  additional  a 
priori  estimates  have  to  ba  obtained  (since  tha  estimates  in  Section  2  only  hold  for 
J  <  0).  for  tha  case  a  3  0  these  estimates  are  given  in  Mock  (1983),  Markovich  (1983). 

Finally,  assume  that  a  is  not  a  functional  on  C*<[0,1])  but  simply  a  function 
(3.8)  a  t  m  ♦  (0#y1#  a  is  continuous  and  nondecreasing,  such  that  the  ionisation  rate 

o(#'(x))  is  space-dependent.  Then  (2.4)-(2.6)  have  to  be  modified  by  substituting  'as* 
r*  M  ■ 

in  the  integrands  by  J  s(+*(s))ds  and  z—  in  the  integration  Intervals  by  that  value 
0  - 

-  x  1 

x  6  [0,1]  for  which  /  a(+*(s))ds  ■  r  holds*  Theorem  3.1  still  holds*  By  estimating 

0 

(2.12)  in  terms  of  Y  an  analogue  of  (2.18),  (2.20)  (also  in  tarms  of  y)  implying 
'♦(1)  ♦  •  as  J  ♦  is  obtained  and  Theorem  3*2  and  Corollary  3*1  follow*  The 

estimates  of  the  current -voltage  characteristic  given  in  Theorem  2*  1  for  0  <  Y  <  J  still 
hold*  The  avalanche-estimate  (2*25)  (with  a  in  the  upper  bound  substituted  by  Y)  holds 
if  for  example  o(t)  >  o^e  with  0^  suff.  large  and  0^  suff.  small* 
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